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Review of Beauty Production at HERA and Elsewhere 
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DESY Hamburg, Germany 



Experimental results on beauty production at HERA are reviewed in the context of 
similar measurements at other colliders. As a result of a phenomenological study of the 
QCD scale dependence of many different NLO and NNLO predictions, a modification of 
the "default" scale choice is advocated. Experimental constraints on the photon-quark 
^ ■ coupling are also investigated. [1] 



1 Introduction 




Beauty production at HERA (Fig. 1) is an important tool to 

investigate our present understanding of the theory of Quantum 

^ ■ Chromo-Dynamics (QCD). On one hand, the large b quark mass, 

D I taken as a hard scale, ensures that the cross sections are always 

Oj' perturbatively calculable. On the other hand, the simultaneous 

D I presence of competing hard scales, such as the transverse momen- 

j-Q^' turn {pt) of the heavy quark, or the virtuality of the exchanged 

photon (Q^), induces additional theoretical uncertainties due to 

^SJ ■ terms in the perturbative expansion which depend logarithmically 

^ ' on the ratio of these scales. The comparison of the measured cross t^- -, ^^ 

.*^^ ■ ^ r igure 1 : Feynman 

sections with theory predictions is therefore particularly sensitive graph for the produc- 
^£^ I to the way the perturbative expansion is made, and can therefore tion of a heavy quark 

potentially discriminate how adequate a particular QCD scheme is pair via the boson-gluon- 
for the dccription of the cross section in question. This can also fusion (BGF) process, 
yield insights for other QCD processes at HERA, and for related 
r~^ I processes at other colliders, including future measurements at the LHC. 

^^ i Since beauty in Deep Inelastic Scattering (DIS) is covered elsewhere |2j, this contribution 

will concentrate on the photoproduction case {Q^ < 1 GcV^), in which the photon is quasi- 
real. For beauty photoproduction at HERA, possible theoretical schemes include 

• The leading order plus parton shower approach, where leading order (LO) QCD matrix 
elements are complemented by parton showers, usually using the DGLAP [3j parton evo- 
lution equations. This approach is implemented in many Monte Carlo models, and mostly 
used for the purpose of acceptance corrections. 

• The kt-factorization approach [1], which can alternatively be used for parton showering, 
combined with the use of generalized parton density functions. 

• The next-to-lcading order (NLO) massive approach [5]. In this approach, the heavy quark 
mass is fully accounted for, and heavy quarks are therefore always produced dynamically 
in the matrix clement, as illustrated by Fig. 1. Alternative LO processes, such as flavour 
excitation in the photon or the proton, are treated as next-to-leading order corrections to 
this BGF process. Processes in which the photon acts as a hadron-like source of light quarks 
or gluons are also included, but make only a small contribution. This approach is expected 
to work best when all relevant hard scales, e.g. px, are of order rrib. 

• For Pt 3> mt, large logpT/nif, terms could in principle spoil the reliability of the pre- 
dictions. In this case, it might be preferable to switch to a so-called massless scheme, in 
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which the b mass is neglected kinematically. The potentiahy large logarithms can then be 
resummed to all orders (next-to leading log or NLL resummation) . Since such an approach 
is obviously not applicable when pT ~ rrib, schemes have been designed which make a conti- 
nous transition between the fixed order (FO) massive, and the NLL massless scheme. This 
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Figure 2: Beauty production cross section measure- 
ments in photoproduction at HERA as function of the 
transverse momentum of the b quark, compared to 
QCD predictions. 



is often referred to as the FONLL scheme [6] 

On the experimental side, several 
different methods arc used to tag the 
beauty final state. The b quark can 
decay semileptonically into a muon or 
electron, which can be identified in the 
detector. The large momentum of the 
lepton transverse to the direction of 
the 6-initiated jet, due to the sizeable 
b mass, can be used to discriminate 
against semileptonic charm decays or 
misidentified light flavour events. The 
finite lifetime of the B hadrons can lead 
to a measureable offset of the decay ver- 
tex with respect to the primary vertex 
of the event, which also leads to a sig- 
nificant impact parameter of the result- 
ing secondary tracks. Finally, a lepton 

tag can e.g. be combined with a lifetime tag, with a second lepton tag, or with a D* meson 
from a b decay. 

Fig. [5] shows a compilation of all recent 
HERA measurements of b photoproduction [7]. 
Reasonable agreement is found with both the 
fixed order NLO QCD prediction [S], and with 
a prediction based on kt factorization [8]. How- 
ever, the data tend to lie somewhat above the 
central prediction in both cases. 

A longstanding apparent discrepancy be- 
tween data and theory in b production at the 
Tevatron was resolved by coinbining a more care- 
ful consideration of B fragmentation and decay 
parameters with an FONLL-based prediction [9] . 
This raises the question whether an FONLL pre- 
dicton, which does not yet exist for b production 
at HERA; would yield an improved agreement. 

For this purpose, consider b production at 
the SppS collider, which had an effective parton- 
parton center-of-mass energy very similar to that 
of HERA. Fig. [3] shows the measured b quark 
cross section [1^ compared to the original NLO 

calculation |11| . Good agreement was observed at a time when the Tevatron experiments 
were starting to claim a discrepancy. Fig. U] shows the same original data [TUI compared to 
the more recent FONLL calculation at b quark and B hadron level [31 [T^] , with identical 
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Figure 3: Integrated beauty production 
cross section at the SppS as function of the 
minimum px of the b quark, compared to 
NLO QCD predictions. 
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parameters as those used for the Tevatron. Good agreement is observed, also at muon level 
[T2] . even though NLO predictions at B hadron and muon level were not available when 
the measurements were made. This indicates that the B fragmentation and decay spectra, 
which had been studied carefully [l0l[T3], were treated consistently in these measurements. 
Furthermore, the NLO and FONLL predictions agree very well with each other, indicating 
that the large logs mentioned above do not yet play a significant role in this pT range (similar 
to the one at HERA). This can also be seen in charm production at HERA [T3] for which 
an FONLL prediction exists. 
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Figure 4: Integrated beauty production cross section at the SppS as function of the minimum pT 
of the B hadron (left) and b quark (right), compared to preliminary FONLL QCD predictions |12] . 

In conclusion, an FONLL prediction for beauty production at HERA would be useful, 
but is not expected to significantly alter the data/theory comparison. 

The dominant contribution to the theoretical error band of Fig. 1 is the variation of the 
renormalization/factorization scale by a factor 2 around the default scale /ig — ■\/m^ + Pr- 
Such a variation is intended to reflect the uncertainty due to uncalculatcd higher orders. It 
might therefore be useful to reconsider this scale choice. 

Ideally, in a QCD calculation to all orders, the result of the perturbatve expansion 
does not depend on the choice of these scalefl In practice, a dependence arises from the 
truncation of the perturbative series. Since this is an artefact of the truncation, rather than 
a physical effect, the optimal scale can not be "measured" from the data. It must thus be 
obtained phenomcnologically. 

Traditionally, there have been several options to choose the "optimal" scale, e.g. 
• The "natural" scale of the process. This is usually taken to be the transverse energy (Et) 
of the jet for jet measurements, the mass m of a heavy particle for the total production cross 
section of this particle, or the combination y^m^ + p^ for differential cross sections of such a 
particle. Often, this is the only option considered. The choice of this natural scale is based 
on common sense, and on the hope that this will minimize the occurrancc of large logs of 
the kind described above, for the central hard process. However, higher order subprocesses 
such as additional gluon radiation often occur at significantly smaller scales, such that this 
choice might not always be optimal. 



''As is common practice, we will not distinguish between the factorization and renormalization scales in 
the following, and set both to be equal. A separate optimization of the two scales, which should bo done in 
principle, will be left for future consideration. 
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• The principle of fastest apparent convergence (FAC) [T5] . The only way to reliably evaluate 
uncalculated higher orders is to actually do the higher order calculation. Unfortunately, this 
is often not possible. Instead, one could hope that a scale choice which makes the leading 
order prediction identical to the next-to-leading order one would also minimize the NNLO 
corrections. Since it can not be proven, this principle, which can be found in many QCD 
textbooks, has not been used very much recently. However, recent actual NNLO calculations 
might indicate that it works phenomenologically after all (see below). 

• The principle of minimal sensitivity (PMS) [16]. The idea is that when the derivative of 
the cross section with respect to the NLO scale variation vanishes, the NNLO crrections will 
presumably also be small. Again, there is no proof that this textbook principle should work, 
but actual NNLO calculations might indicate that it does (see below). 

To illustrate these principles, consider two examples. First, the prediction for the total 
cross section for beauty production at HERA-B [17] (Fig. [5]). The natural scale for this 
case is the b quark mass, fiQ = mi,, and all scales are expressed as a fraction of this reference 
scale. Inspecting Fig. [S] one finds that both the PMS and FAC principles, applied to the 
NLO prediction and to the comparison with LO (NLO stability), would yield an optimal 
scale of about half the natural scale. The same conclusion would be obtained by using the 
NLO-fNLL prediction, including resummation, and comparing it to either the LO or the 
NLO prediction (NLO+NLL stability). 
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Figure 5: Scale dependence of the total cross section for beauty production at HERA-B [17] (left) 
and for Higgs production at LHC [18| for two different masses (right). 

Second, the prediction for Higgs production at the LHC [TS] (Fig. \5^. The reference 
scale is now the Higgs mass (/ig = ttih)- However, inspecting the behaviour of the LO 
and NLO predictions, neither the FAC nor the PMS principle would yield a useful result 
in this case, since the two predictions do not cross, and the NLO prediction does not have 
a maximaum or minimum. This situation occurs rather frequently, and is also true for b 
production at HERA. Fortunately, in the case of Higgs production, the NNLO and even 
NNNLO predictions have actually been calculated (Fig [5|) . Applying the FAC and PMS 
prescriptions to these instead (NNLO stability), again a scale significantly lower than the 
default scale would be favoured. This might indicate that choosing a scale which is sinaller 
than the default one makes sense even if the FAC and PMS principles do not yield useful 
values at NLO. 
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Beyond these examples, a more general study is needed to phenomenologically validate 
this approach. To avoid additional complications arising from a multiple scale problem 
caused by e.g. the scale Q^ at HERA or the scale AIz at LEP, the study was limited to 
cross sections for photoproduction at HERA, or hadroproduction at fixed target energies, 
the Tevatron, and LHC. The somewhat arbitrary selection of processes includes beauty 
production at the SppS [l9l[20|, the Tevatron [20], and HERA-B [17], top production at 
the Tevatron [T71 HD] , direct photon production at fixed target [H] , Z [52] and Higgs [TS] 
production at the LHC, jets at HERA [23] and at the Tevatron [24]. This selection is 
obviously not complete. However, it is not biased in the sense that all processes that were 
considered were included, and none were discarded. 

In each case the natural scale as defined 
above was used as a reference. In addition, 
wherever possible, the optimal scales from both 
the FAC and PMS principles, evaluated at NLO 
(NLO stabiUty), NLO+NLL (NLO+NLL stabil- 
ity), and/or NNLO/NNNLO (NNLO stability) 
were evaluated separately. Fig. [Hj shows the 
result of this evaluation. Each crossing point, 
maximum, or minimum in Fig. [5| yields one en- 
try into this figure, and similarly for all the other 
processes. The conclusion is that the FAC and 
PMS principles tend to favour scales which arc 
around 25-60% of the natural scale. Amazingly, 
this seems to be independent of whether these 
principles are applied at NLO, NLO+NLL, or 
NNLO level. For the jet |24] or b-jet [25] cross 
sections at the Tevatron, it has in part already 
become customary to use half the natural scale 
as the central scale. 

Using the natural scale as the default and varying it by a factor two, which is the choice 
adopted for most data/theory comparisons, covers only about half the entries, while the 
other half lies entirely below this range. Instead, using half the natural scale as the default 
and varying it by a factor two, thus still including the natural scale in the variation, covers 
about 95% of all the entries. 

This yields the following conclusions. 

• Obviously, whenever an NNLO calculation is available, it should be used. 

• Whenever possible, a dedicated scale study should be made for each process for the kine- 
matic range in question. Although there is no proof that the FAC and PMS principles 
should work, in practice they seem to give self-consistent and almost universal answers for 
processes at fixed target energies, HERA, the Tevatron, and the LHC. 

• In the absence of cither of the above, the default scale should be chosen to be hali the 
natural scale, rather than the natural scale, in particular before claiming a discrepancy 
between data and theory. Empirically, this should enhance the chance that the NNLO 
calculation, when it becomes available, will actually lie within the quoted error band. 



_L 



0,4 0,5 0,e 0,7 0,81 



scale ^7^ 



Figure 6: Summary of optimized scales de- 
rived as described in the text. 



DIS2007 



. 1 1 . 1 . 1 1 1 1 1 1 . 1 

^ , t.^ . 


1 . 1 . 1 1 11 1 1 1 1 1 1 1 

do/dp^(ep^ebX) "i 
Q^<1GeV^ 0.2<y<0.8, hV 2 


B % 


[ 1 NLOQCD(FMNR) 


■ I * H1 99-00 b jet ■■■J.,it| 
"/ T H1 99-00 b^n let ttf. 
f • ZEUS 96-00 b^n jet 

; A ZEUS (prel) 04 b->n jet 

■ O ZEUS 96-97 b^e 


■■■■■■1. 1 
'■■■■*-..,* 

■■■f......... 


r D H1 97-00 b^D*n 


■ ZEUS 96-00 b^D'n 

* ZEUS (prel) 96-00 bb^jin 

, ... 1 .... 1 ... . 


>■■■: 
1 , ... 1 , ... 1 ... . 



HERA 

Now consider the application of the ntn^A 

last proposal to actual data/theory ^10^ 

comparisons. Fig. [7]shows the resulting § 

comparison for beauty production at ~ ^ 

HERA. Although before the change of it 

the default scale the agreement was al- "B 

ready quite reasonable, the new choice, ^° 

based on theoretical/phenomenological 

arguments, improves the agreement. A 1 

similar statement [T] can qualitatively 

be made for beauty production at the ^p' 

SppS [TU], the Tevatron [55], and even _, b^ iGeV) 

at RHIC [27], where a discrepancy has ^. _ 

, 1 • 1 mi , -,1 rigurc 7: Beauty production cross section measure- 

been claimed. ihe agreement with °^ . , ^ , i- ^ TTT^n a r ^- r^i 

"Ij i__j. ments m photoproduction at HERA as function or the 

^ ' — ! transverse momentum of the b quark, compared to 

as charm [28] and top [29] production qqd predictions 

at the Tevatron also improves, or at 

least does not get worse. The same is true for inclusive jets at HERA, both in DIS [50] 

and in photoproduction [23]. In one case [31] half the scale has already been used in a 

published HERA result. 

So far, no example is known to the author where the proposed change of default scale 
would result in a significant worsening of the data/theory agreement in a photo- or hadropro- 
duction cross section. Thus, the phenomenologically motivated change seems to be sup- 
ported by the data. It should therefore be considered to make it the default for future 
predictions at HERA, the Tevatron, and the LHC. 

The investigations of the photon quark coupling are reported elsewhere [1] . 
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